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Abstract: We consider a supersymmetric SO(IO) model with a SU(3) symmetry of flavour 
in which fermion masses emerge via the see-saw mixing with superheavy fermions in 16 + 16 
representations. In this model the dangerous D = 5 operators of proton decay are natu- 
rally suppressed and flavour-changing supersymmetric effects are under control. The mass 
matrices for all fermion types (up and down quarks, charged leptons as well as neutri- 
nos) appear in the form of combinations of three rank-1 matrices, common to all types of 
fermions, with different coefficients that are successive powers of small parameters, related 
to each other by SO(IO) symmetry properties. Two versions of the model are considered, 
in which approximate grand unification of masses takes place between quarks and leptons 
of the first family (with very small tan /3) or for the ones of the second family (predicting 
moderate tan/3 ~ 7-8). The second version exhibits an interesting mechanism of unifica- 
tion of the determinants of the Yukawa matrices of all types of fermions at the GUT scale 
and it provides a perfect fit of the known data for fermion masses, mixing and CP-violation. 
It predicts a hierarchical pattern of neutrino masses with non-zero 9e3, within 2-7 degrees. 
In addition, it predicts the correct sign of the baryon asymmetry of the Universe via the 
leptogenesys scenario. 
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1. Introduction 

In the recent years the experimental information on fermion masses and mixing has become 
increasingly accurate, with hierarchies and apparent regularities representing a puzzle that 
continue to call for new theoretical frameworks beyond the Standard Model. Among these, 
the most promising possibility is related to the minimal supersymmetric extension of the 
Standard Model (MSSM). In the MSSM, the fermion masses and mixings emerge from the 
Yukawa couplings to the Higgs fields h^: 

qXuu" hu + qXdd'' + IXee"" + IXy^v" , (1.1) 

where q = (n, d), I = {u, e) are weak isodoublets and u'^, (F, e^, v'^ are isosinglets (the family 
indices are suppressed) and A's are matrices in family space. ^ Also the neutrino Yukawa 

^Throughout the paper we will represent quantities that are 3x3 matrices of flavour with a hat, i.e. A — 
{Xij}, and denote flavour indices with i, j. 
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couplings are included that induce their Dirac mass terms. However since v'^ are neutral, 
they can have Majorana mass terms M i^'^Xycv'^ where M is a. large mass scale related to 
some physics beyond the Standard Model. As a result the effective D=5 operators are 
induced by the seesaw mechanism [j^ : 

-j-^ II huhu , Xy = XujyXjy^Xl,^ . (1-2) 

Regarding charged fermions, we observe that their Yukawa couplings (eigenvalues of 
the matrices A^^^^e) show strong hierarchy in family space: 
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: Xc 


: Ai 
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: en : 1 


en ~ 1/500 
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■ ed : I 
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: ee : 1 


ee ~ 1/50 


(1.3) 



while the quarks mixing angles are small: si2 ~ S23 ~ ed, S13 ~ e^. In the "up" sector 
(U) the scaling pattern ( |1.3| ) is almost exact while in the "down" and "charged leptons" 
sectors (D, E) there are sensible deviations. The deviations are present mainly for the first 
two generations in the D, E sectors, and we will naturally connect this with the largeness 
of the Cabibbo mixing angle, as compared to the other quark mixings. Prom the above list 
we note that the "down" and "charged leptons" sectors have similar hierarchies, reaching 
the well known observation that e^ ^ eg — e^. 

On the other hand, in the neutrino sector, the neutrino mass eigenstates show a milder 
hierarchy (they could also be strongly degenerate) while the mixing angles are large. 

The concept of grand unified models, together with the idea of family unification pro- 
vides a promising framework for solving the flavour problems. In this respect, probably 
the SO(IO) is the most interesting candidate, since it unifies all the fermions of one family 
into a single multiplet 16 and so it can naturally link the Yukawa constants of different 
fermion sectors, with typical Clebsch factors of SO(IO) |^. For family unification the max- 
imal symmetry group is SU(3) (or U(3)) which unifies the three families in one horizontal 
triplet [P, In this way, the inter-family hierarchy can be related to the SU(3) breaking 
pattern. 

In the context of SO(IO) xSU(3), the three families of fermions fit into the representa- 
tion (16,3) = 16*, where i = 1,2,3 is a family SU(3) index. The family symmetry forbids 
them to have Yukawa couplings with the Higgs 10-plet of SO(IO), since this is singlet of 
flavour while the bilinear of fermions transforms as 3 x 3 = 3 -|- 6. Therefore, in order to 
generate fermion masses one needs a higher order operator of the form Gjj 16*16-^10, where 
Gij is a field-dependent effective Yukawa coupling, that transforms as 3 x 3 under SU(3), 
for example it may be an (anti)sextet. It is clear that the pattern of fermion mass matri- 
ces will reflect the VEV structure of G i.e. the horizontal symmetry breaking. However, 
this operator should also be SO(IO) dependent, otherwise it would leave in each family an 
undesirable degeneracy between up, down and charged lepton masses. Therefore one is led 
to consider operators of the form 

Gij{X)16'16nO , (1.4) 
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where X is a set of SO(IO) fields, and we will see that one can efficiently take 45-plets 
of SO (10). These higher order operators can be effectively generated from renormalizable 
couplings involving heavy vector-like fermions, so called universal-seesaw mechanism |^]. 

In this paper, we follow the idea |^, |^ that the Yukawa matrices are built as a linear 
combination of three fixed rank-one projectors in flavour space, i.e. three 3x3 symmetric 
rank-one matrices A, 2, 3:^ 

Xf (xejPi + efp2+ P3, f = u,d,e,i'D ■ (1-5) 

Only the coefficients depend on the fermion type (/ = u, d, e, vo) while the three disoriented 
rank-one projectors are common to all types. 

We will build models where the matrices Pn are generated from the VEVs of flavon 
scalar fields breaking the SU(3) symmetry, and the factors e/ arise from the VEV of the 
X fields at SO(IO) breaking and are related by specific Clebsch relations. 

As far as the Pn are rank one matrices with generically order one entries, the eigenvalues 
will roughly follow the hierarchies (^]^). We will see that these coefficients are realized 
naturally using the VEV of a single 45x field with a generic orientation compatible with 
the SO(IO) breaking to the Standard Model group. Then the breaking of flavour realized by 
the Pn projectors not only generates the flavour mixing angles, but also the right deviations 
of the eigenvalues from exact hierarchy. Namely, there will be a link between the deviation 
in the D, E sectors and the largeness of the Cabibbo angle. 

We will provide the model realizations of this ansatz that also satisfy the naturality 
conditions, namely the dangerous D=5 operators inducing proton decay are automati- 
cally suppressed and the supersymmetric flavour changing effects are under control. These 
models will lead to three different cases of the general ansatz (|1.5| ) for the fermion mass 
matrices. One of these cases will be found in perfect agreement with the observed masses, 
mixing angles and CP violation, and quite interestingly this case points to an exact grand 
unification of determinants of the Yukawa matrices at the GUT scale: 

det \u = det = det Ae = det X^j^ ■ (1-6) 

The paper is organized as follows: in section ^ we introduce the "universal seesaw" 
mechanism for fermion masses, and discuss the VEV patterns for the SO(IO) and SU(3) 
symmetry breaking. In section ^ we describe the model realizations in terms of renormaliz- 
able operators and discover three interesting cases corresponding to the Yukawa couplings 
unification in the first, second or third fermion family. Singling out the second case, in 
section |5| we will describe the unification of Yukawa determinants and its predictions for 
fermion masses and the tan fj parameter. In section ^ we will solve analytically the model 
in terms of its parameters and also perform a numeric fit. 

A preliminary short version of this work has been presented in HI. 

^One may choose also non symmetric rank-one projectors, but at the price of loosing predictivity. 
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2. Universal seesaw in SO(IO) and flavour symmetry 

Along the lines of Q , we describe the generation of effective Yukawa couplings at the GUT 
scale via the universal seesaw mechanism that employs SO(IO) gauge symmetry together 
with the SU(3) horizontal symmetry. 

The model uses Higgs fields in the following representations of SO(IO): 16, 16, 10, 
45, 54.^ For later reference we recall their decomposition in terms of the standard SU(5) 
and Pati-Salam SU(4) xSU(2) xSU(2) subgroups: 

16 = 1 + 5 + 10 = (4, 2, 1) + (4, 1, 2) , 
T6 = l + 5+T0 = (4,2,l) + (4,l,2), 
10 = 5 + 5 = (6,1,1) + (1,2,2), 

45 = 1 + 10 + To + 24 = (15, 1, 1) + (1, 1, 3) + (1, 3, 1) + (6, 2, 2) , 
54 = 15 + 15 + 24 = (1, 1, 1) + (1, 3, 3) + (20', 1, 1) + (6, 2, 2). 

We assume that SO(IO) is broken to the Standard Model group in a single step by 
means of the set of Higgs fields 16h + 16h, 54/^ and 45, with VEVs of the same order. In 
particular the 16h and 16//, with VEV towards lsu(5), breaks SO(IO) down to SU(5), while 
the 54 with VEV oc diag(l, 1, 1, -3/2, -3/2) breaks SO(IO) down to SU(4)xSU(2)xSU(2). 
The intersection of these two breaking channels leads to the Standard Model symmetry 
group SU(3)xSU(2)xU(l). In addition, internal consistency of the Higgs sector requires 
the presence of fields in the 45 representation. 

We assume that we have three 45-plets, 45^/,, 45r, 45x- The first, with direction 
towards the (15, 1, 1) component, {45bl) oc (1, 1, 1,0,0) a, can be used for the solution 
of the doublet-triplet splitting problem via the Missing VEV mechanism ^j. A second 
one, with VEV orthogonal to that of 45bl, towards the (1,1,3) component: (45j?) oc 
(0, 0, 0, 1, 1) a. Finally we have a 45x with generic VEV towards both these directions.'^ 

The Higgs doublets {hu,hd) = </>, which induces the electroweak breaking and the 
fermion masses is sitting in the (1, 2, 2) component of the 10-plet. The (6, 1, 1) component 
contains the colored higgses T, T, that should be very heavy, of the order of the GUT 
scale, in order to suppress proton decay via dimension 5 operators p^ . 

As far as the breaking of the flavour SU(3) symmetry is concerned, we will use three 
sextets Xi2 3 with disoriented rank-one VEVs H, H]. These VEVs generically can be 
presented as (xn) = Un diag(0, 0, 1) C/* , with Un being 3x3 unitary matrices. Alternatively 
a sextet field Xn may be built as effective tensor product of triplets x*"* = C ' C'^ [Q- If 
one introduces three triplets (n = 1,2,3) having misaligned VEVs in SU(3) space, then 
the tensor products • ^.^ automatically provide effective Xn with VEVs of rank-one. 

^Motivated by minimality arguments, we avoid using huge representations like 210, 126 + 126 etc. . . 

*The superpotential(s) structure realizing the VEV pattern and a detailed analysis can be found in p^ , 
PH , p^ , and the needed form of the potential may be motivated with an additional discrete symmetry or 
just rely on the nonrenormalization theorem. The essence of the Missing VEV mechanism also motivates 
the use of the three adjoint fields 45 x, 45 r, 45bl- 
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Figure 1: Universal seesaw mechanism for Yukawa couplings. 



2.1 Universal seesaw 

We introduce heavy vector-like fermions 16'-, 16[ that are antitriplets of flavour and consider 
the following superpotential terms, allowed by the SO(10)xSU(3) symmetry: 

1^ = 71016*16^ + ys 16*16; +T6:M^^^^^yl6;. (2.1) 

where 10 is the fundamental Higgs and S is a Higgs field that could be in the 1, 45 or 210 
representations and whose role will be exploited in the next section. The flavour structure 
is encoded in Mheavyj while the first two terms are flavour universal. Also, Mheavy is still 
an effective SO(IO) operator. 

Near the GUT scale where some higgses develop a VEV, the heavy fermions 16' get 
mass and decouple, so the light fermions 16 acquire an effective Yukawa-like coupling with 
10 that is approximately given by the SO(IO) "seesaw" formula: 

l^Hght = G,, 10 16' 16^' , G = 7 y (S) M-!^,y . (2.2) 

We begin to see from here that the Yukawa couplings of the light fermions will be propor- 
tional to the inverse of the heavy fermions mass matrices. 

To see this happen in detail, it is useful to decompose the fermion fields under the 
Pati-Salam group as 16* = /*(4,2,1) + /^'(4,1,2) and 16^ = J^i(4,2,l) + i^^(4,l,2), 
T6. = J"f (4, 2, 1) + Fi(4, 1, 2). Note that / = g, / and = q', I' are weak isodoublets while 
= u'^, (F^ e'^, u^, F = u' , d' , e', z^'., are weak isosinglets. With this decomposition, we can 
illustrate the couplings present in (2.1), before seesawing, as: 

jc pc j:c 



f 
F 

J" 



70 /Ujr ^ 

HF Mp 
\j(t> M^J 



(2.3) 



where Mheavy is decomposed in Mp and Afjr , respectively the mass matrices for (unmixed) 
isosinglets and isodoublets (recall that heavy fermions F-F'^ and T-T'^ have vector-like 
masses) and /ii?, /ijr are the projections of 7'(S) on the isosinglets and isodoublets chan- 
nels. Recall that the latter entries are flavour blind and the nontrivial flavour structure is 
contained only in the matrices Mp and Mj^. 
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Figure 2: Universal seesaw mechanism for RH neutrino mass.. 

This form makes explicit the seesaw mechanism that follows from the mixing between 
the isosinglets and F'^ and between the isodoublets / and J-. After this mixing one 
ends up with an effective Yukawa couplings of the light fermions /, to the Higgs 0, as 
in the Standard Model: (p f Xfcf^. Decomposing further the SU(4) of Pati-Salam in quark 
and lepton channels one recovers the Standard-Model Yukawa couplings with Yukawa 
matrices given by:^ 

Xf = jfiFMp^ +lf^rM^\ f = u,d,e,UD. (2.4) 

We see that the Yukawa couplings receive two contributions, from the inverses of both Mp^ 
and M^^. Therefore each fermion type u, d, u, e will receive one contribution from the 
mass matrix of its relative heavy isosinglets F = u', d', v' , e', and an other from that of 
the appropriate heavy isodoublets F = q' I'. We will expand below the form of Mp and 
Mjp to describe explicitly the realization of the charged fermions mass matrices. Before 
doing this, we complete this section addressing the generation of neutrino masses. 

The Yukawa couplings obtained above generate Dirac masses for charged fermions as 
well as for neutrinos and, as is frequent in unified models, these neutrino masses will be 
unrealistically high. Therefore one is led to consider also Majorana mass terms for the 
RH neutrinos, and to generate the LH neutrino masses by canonical seesaw. After all RH 
neutrino are neutral particles and there is no reason to forbid their Majorana mass term. 

The Majorana mass for the RH neutrinos can be generated in a similar way as above, 
via an universal seesaw mechanism using the 16h Higgs field, whose VEV selects the RH 
neutrino i^ji from the 16-plets. The relevant part of the superpotential is in this case: 

Wn = 7nT6h 16' 1, + M^' hlj (2.5) 

where 16h is the Higgs whose VEV preserves SU(5) as described above, and where we have 
introduced a flavour triplet of SO(IO) singlet fermions (1, 3) = 1,.^ The flavour structure is 
again encoded in the mass matrix of the heavy singlets, while the first term is flavour 
universal. Also, M^v is again an effective operator that should of course be an SO(IO) 
singlet. We do not want at this point to dwell on the details of its construction, that will 
follow lines similar to the case of the charged fermions, described in the next section. 

'""Without the risk of confusing notations we denote here with / the fermion type, instead of the hght 
isodoublets. 

^Alternatively one could use e.g. a (45, 3) multiplet. 
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When the heavy singlets Ij get mass and the 16h develops a VEV, the RH neutrino 
mass matrix emerges from the universal seesaw, as depicted in figure ^, and is given by 

M,^ = C^M^^ , (2.6) 

where C = 77v(16//). Then, the light (LH) neutrino mass matrix results from the canonical 



seesaw with the neutrino Dirac couplings given by (2^), as 



M,, = M,Jl-^Mi^ = l^M.^M^Ml^ . (2.7) 

with M^^ = VuXuD where Vu is the "up" electroweak VEV, Vu = {hu) = f sin/3 ~ 
174 GeV sin/3 in the MSSM. 

We can at this point comment on the scales required for this neutrino mass generation, 
that can be inferred from relation ( |2.7D . Since the Dirac neutrino masses are unified with 
the other fermions, they will fall in the GeV range; in particular in these models Mjy^ ~ 
0.01-10 GeV (see later). Therefore in order to obtain realistic neutrino masses of magnitude 
<0.1eV one requires the two scales C, Mn to be related as (C/GeV)^ ~ 10^2(MAr/GeV). 
This leads for example, for Mn ~ 10^^ GeV, to C ~ 10^^ GeV. Since C = 77v(16h), this 
requires ^jy to be just of order 10^^, in order to keep (16^^) at the GUT scale. On the 
other hand there is no problem in having such a high M^v decoupling scale, since the Ij 
fermions are SO(IO) singlets and they do not interfere with the gauge group breaking. 

Finally, for later convenience, we observe that introducing dimensionless matrices also 
for the left and right handed neutrinos )^ul,uii = M^^^^^/vu and for the heavy singlets 
Xn = Vu/C^, we can write the above expressions in the form: 

2.2 Flavour structure from rank-one decomposition 

Let us discuss now the form of Mheavy Under flavour SU(3) symmetry, Mheavy transforms 
as 3 X 3, therefore it can be generated via the VEVs of some flavon flelds in SU(3) sextet 
representation. 

Our central assumption now is that Mheavy is built as a linear combination, with 
hierarchic coefficients, of three rank-one matrices in flavour space. The sum will give a non 
degenerate matrix, with hierarchic eigenvalues. 

The three rank-one matrices may for example be built out of sextet or triplets scalar 
flelds. For example there may be three sextets, that break the flavour symmetry by devel- 
oping a VEV, and their VEVs should be of rank-one. This can be achieved through some 
mechanism like those described in [15|. 



The hierarchic coefflcients will be generated as effective operators at SO(IO) breaking, 
and we will focus on them later; at this point it is useful to denote them with 0!]jeavy,n! 
with n = 1,2,3: they are singlets under flavour but are functions of SO(IO) fields X and 
generically transform as 1, 45 or 210. 

Summing up, we introduce three sextets x*"* ^-nd mix them with the Q^heavy's* 

^h^avy = E «heavy,„(X)x^ (2.9) 
n=l,2,3 
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Even the heavy singlets mass matrix that gives Majorana mass to neutrinos transforms 
as a sextet of flavour, and thus it may be naturally built using the same sextets Xn , however 
with different coefficients ajyn- 



E 

?i=l,2,3 



(2.10) 



At the GUT scale, where SO(IO) is broken to the SM group, the coefficients Oheavy,n 
decompose in the weak isosinglets and isodoublets channels, and consequently Mheavy de- 
composes in Mp and Mjr. We can write, for all the heavy fermions 

Mp = M[aF,i Qi + aF,2 Q2 + of.s Q3] , 
Mjr = M[ajr^i Qi + aj^ i Q2 + a^^i Q3] , 

Mn = M[aN,i Qi + aN,2 Q2 + aiv,3 Qs] , (2.11) 

where the a's have been projected on the isosinglets or isodoublets sectors and are now 
pure numbers, and where we have parametrized the rank-one VEV of the sextets (xn) with 
MQn- M is the scale of flavour breaking and Qn are dimensionless flavour projectors. 

To find the Yukawa couplings and the RH neutrino masses, given by the inverses of 



these heavy mass matrices via (2.4) and (2.6), we can exploit a useful property of every 
combination of rank-one projectors, namely that its inverse is also a combination, with 
inverse coefficients, of (new) rank-one projectors. Explicitly, 



(aiQi + a2Q2 + a^Qs) ' = ( —A + — A + — A 

ai 02 03 



(2.12) 



where the Pn are the "reciprocal" of the Qn projectors: if one parametrizes the Qn with the 
symmetric product of generic complex vectors g„ (flavour triplets) as Qn = QnQn = Qn'^Qm 
the Pn can parametrized with three vectors p„, that are in fact the reciprocal of the g„ 
ones: Pn = ^^nmrQm A Qr/Dq, with Dq = Qi ■ (^2 A Qs). In this notation, the inverse heavy 
mass matrices are then: 

= M-i [a-\ A + A + A] , 
M^i = M-^[a^\ A + A + A] , 



M 



N 



[a]^\ A + A + a^ig A] 



(2.13) 



From (|2.4| ), (|2.6f ) we then find the form of Yukawa couplings and RH neutrino masses, 
reaching the conclusion that they are all written as combinations of the same three rank-one 
flavour projectors Pn- 



Xf oc 



Kr oc 



Pi + a/,2 P2 + a/,3 -^3 



P2 + ^^3 



f = u,d, e, UD 



(2.14) 
(2.15) 
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with coefficients tfiat are 



a 



VR,n 



aF,n 

0:N,n 



(2.16) 
(2.17) 



As far as tlie liglit (LH) neutrino Major ana mass matrix is concerned, we finally have 
to perform the canonical seesaw ( |2.8| ), and we can exploit an other interesting property of 
rank-one combinations. Indeed, the following "seesaw" relation holds: 



(ai Pi + 02 P2 + 03 P3)(/5i Pi + P2P2 + Ps P3) (ai Pi + 02 P2 + as P3 



«i p , «2 p , "3 p 

f3i (32 



(2.18) 



i.e. the seesaw acts on the coefficients only, but does not change the flavour projectors 
Pn- Hence, also the neutrino mass matrix is built as a combination of the same flavour 
projectors: 

^,2 



A 



OiuL,l Pi + auL,2 P2 + P3 



a 



a 



I'D 



a 



(2.19) 



We have thus shown that the basic idea of dealing with combinations of common rank- 
one projectors is not spoiled by the seesaw generation of neutrino masses, and this will 
allow below the model to be nicely predictive in this sector. 

2.3 Choice of 45r in the "right" direction 

In the universal seesaw for Yukawa couplings, one can take the Higgs field S in the 45 
representation. Then, the possible VEV directions are (15,1,1) and (1,1,3) of Pati-Salam. 
Choosing the (1,1,3) direction leads to a number of interesting consequences, that we 
will briefly describe. Let's denote S in this case as 45/j. When its VEV is chosen in that 
direction, i.e. 7'(45j?) = /i_F(l, Ij 3), we have that it couples only light and heavy isosinglets, 
but not isodoublets. Correspondingly iijr = and in the couplings matrix ( |2.3| ) the f-J^'^ 
entry is zero: 

jc pc -pc 

/ / 70 \ 
T\jc^ Mr J 



(2.20) 



The first consequence to be noted is that now only the f^, participate in the seesaw, 
and instead of ( ^.41) , it yields a simpler expression for the Yukawa matrices: 



(2.21) 



This expression has the property that a given hierarchy of the heavy isosinglets fermions 
F is directly reproduced in the light Yukawa couplings, while in the general case the light 
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Yukawa was given by the sum of the inverses Mp^, M^^. Each "hght" fermion type / 
receives now a Yukawa coupling only from the "heavy" relative sector of isosinglets F, 
proportional to their inverse mass matrix in flavour space. This is a particularly desirable 
situation, since every fermion sector follows an almost exactly hierarchic pattern, that is 
easier to obtain from a single inverse matrix, rather than the sum of two. In this sce- 
nario therefore, we observe that the heavy F fermions should exhibit an inverted hierarchy 
pattern, with the lightest particle being the heavy correspondent of the "top" . 

In addition to these facts, we note also two other remarkable features of the universal 
seesaw using a 45ij: 1) the LLLL (dominant) part of the D=5 proton decay is automatically 



eliminated, because (1,1,3) does not couple the light and heavy isodoublets [g, |10[; 2) the 
sfermion mass matrices are automatically aligned with the square of the fermion ones, thus 
avoiding the SUSY flavour problem (see e.g. [0). We conclude that the breaking in the 
SU(2)jj direction automatically allows the correct mass generation, suppresses the proton 
decay and avoids the SUSY flavour problems. 

2.4 Flavour Clebsches from 45x 

We focus now on the SO(IO) structure for the construction of the hierarchic coefficients 
«F,i,2,3- We will argue that the use of a 45x Higgs field alone predicts hierarchy parameters 
that are in agreement with the phenomenologically observed pattern. 

We first observe that since the known fermion masses follow a direct hierarchy, one 
needs for the heavy fermions an inverted pattern, i.e. a_F,i,2,3 ~ liCie^- We may generate 
small parameters from the ratio of a VEV (45x) to some higher scale V , and the sequence 
can be generated by taking successive powers. For example the operator form of M^eavy 
may be: 



Keavy ^ 



(2.22) 



whose realization as effective operator will be addressed in detail in the next section. 

At SO(IO) breaking, when 45x develops a VEV, it is clear that the three coefficients 
will be hierarchic according to the small ratio (45j)s:)/^, and this, projected on the different 
components of the heavy 16-plets, will give rise to different small parameters. Let's denote 
them ed, Ce, ei/^ for isosinglets and e^, ei for isodoublets. When Mheavy is decomposed 
in Mp^ Mjr and all fields get their VEV, we have: 

Mi. oc [Qi + eF<02 + , F = u,d,e,UD (2.23) 

Mr « [Qi + e^Qs + e^Qs] , J^ = q,l. (2.24) 

The Yukawa couplings will be proportional to the inverse of the isosinglets matrix Mp 



(see ( 2.21 )) leading to: 

A/ oc 



A + e/ A + A , f = u,d,e,i^D, (2.25) 



that inherits the small hierarchy parameters from Mp, but in inverted order. 
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The ef parameters are determined by the breaking pattern of 45 x, so an important 
question is: can they be reaUstic? The answer to this question is not only that they can, 
but also that this fact is a prediction when one uses the 45 x- To see this, let us parametrize 
the generic allowed direction of {45x) as {45 x) /V = ei5(15, 1,1) +€3(1,1,8). Then the 
different small parameters e^, eg, e^^^, eg, ei are all determined by the two eis, £3 with 
some SO(IO) Clebsch-Gordan factors, as follows: 

= ei5 + £3 = — 3ei5 + €3 

erf = ei5 - es = -3ei5 - £3 (2.26) 

£(? = — €15 = 3ei5 . 

Hence, since they depend on just two parameters, they are not independent, but satisfy 
the following relations: 

3 

ee = -ed-2en, e,y^ = -2erf - e„ , = -3eg = -(e^ + e„) . (2.27) 

The first relation is particularly interesting, since it leads to the following two obser- 
vations. Firstly, once one ensures that |e„| <C je^l, the equality of hierarchies in the D, E 
sectors is predicted: |ee| — |ed|. Therefore the approximate D-E symmetry appears as a 
prediction, assuming just that the U sector is more hierarchic. Secondly, still from the first 
relation, we observe that the signs of eg, are nearly opposite, and this will be crucial 
for explaining the deviations from exact hierarchies in the E and D sectors, in section |6|. 
Therefore we conclude that that the use of a 45x with VEV in a generic direction com- 
patible with SO(IO) breaking to SM naturally accommodates the pattern of hierarchies 
observed in nature. 

Regarding the neutrino sector, one similarly introduces an operator for the singlet 
states built with the same flavon sextets, that when inverted generates a RH neutrino 
mass matrix in the form of a combination of the same Pn projectors, with a different 
hierarchic parameter e^j^. 

The only important difference with respect to the case of Yukawa couplings is that 
the three coupling constants the were hidden in the previous expressions (like for exam- 



ple ( 2.22 )) now play a role: indeed there one was able to reabsorb the coupling constants 
in the normalization of the Xn sextet fields, but this can be done only once. Now one can 
only absorb one in the overall scale and an other in the hierarchic parameter e^^ , while the 
third will remain as an explicit parameter, a: 

1 



a 



(2.28) 



This additional constant will appear also in the mass matrix of the light neutrinos: in 
fact using X^^ in the seesaw with the Dirac neutrino matrix ( |2.25| ) and exploiting for the 
rank-1 flavour decomposition the seesaw relation ( 2.18 ), ( 2.191) , we find the form of the LH 
neutrino mass matrix: 



el^ Pi + e^^ P2 + aPs 



£2 

with eu^ = (2.29) 
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We have found that the light neutrino mass matrix finahy has the same form of combi- 
nation of the three Pn rank-1 projectors, but with non-hierarchic coefficients. The complex 
numbers 6,^^ and a will be the two free parameters to be used when testing the model 
predictions in the neutrino sector. 



3. Three Yukawa unification cases 



Coming back to Yukawa matrices, the combination ( ^.25 ) has a/,1,2,3 = ^ f • ^/ • 1) o'^s 



can observe that, neglecting the corrections of order Cf, the fermions of the third generation 
will roughly have all the same Yukawa coupling, corresponding to the eigenvalue of P3 . This 
realizes the traditional t-b-r unification scheme, where at GUT scale the Yukawa of the 
third generation get unified. In the MSSM this requires the parameter tan /? = Vu/va to 
have a very large value, ~ 55. 

However, we note that the pattern of known fermion masses does not require this to 
happen, it just requires the three coefficients to be hierarchic within each fermion sector. 
Therefore one may add an overall power of to the construction of the last section and 
look for Yukawa couplings of the form: 



ejPi +efp2 +P3 



f = u,d,e,UD (3.1) 



where the power of ej that we have added is a generic integer because it should be generated 
as effective powers of 45x/V- 

This form, for n = 0,1,2, singles out three interesting cases of approximate Yukawa 
unifications, that corresponds in the MSSM to the cases illustrated in figure ^, for different 
values of tan /?: 

I) for n = 0, the coefficients are 1, ej^ , ej'^, and the GUT Yukawa couplings of the first 
generation are approximately unified. In the MSSM this requires a very low value of 
tan/? (< 2); 

II) for n = 1, the coefficients are ef,l,ej^, we find an approximate unification of the 
second generation, at moderate tan/? (~ 10); 

III) for n = 2, the coefficients are e^, ej, 1, and large tan f3 (> 50), one reaches unification 
of the third generation, as in the known case of t-b-r unification. 

In figure ^ we have plotted the Yukawa couplings calculated from experimentally known 
masses, assuming MSSM renormalization, and including experimental uncertainties as well 
as uncertainty from large At renormalization. 

One may observe that in the first two cases unification of the first or second generation 
Yukawa couplings is not exact (there is a splitting of a factor of 4 in the ratios Xd/^e or 
As/A^); This is due to yukawa eigenvalues deviating from exact hierarchy (a straight line) 
and we will see that this deviation is explained, in the present framework, as a consequence 
of the magnitude of the Cabibbo mixing. We will also show that among the three cases the 
second one (II) gives optimal agreement with the data. In this case, also, the intermediate 
tan/3 is in the right range, favored by the MSSM Higgs sector constraints |18]. 
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tan 13^2 tan /3^7 tan ^^50 




Figure 3: Yukawa couplings at GUT scale calculated from experimental charged fermions masses, 
assuming MSSM running and including uncertainty from experiments and from At renormalization. 
U, D, E are respectively red (continuous), green (dashed) and blue (dot-dashed). 



4. SO(IO) realizations 



To realize these three models of Yukawa couplings, for generic n, one should look for 
effective heavy fermion mass operators of the form. 



45 



X 



Xi 



, 45x ij 

+ -i—X2 



+ 



452^ 



in 

2 ^3 



(4.1) 



where again the x*/2 3 sextets of flavour with VEV that gives the <3i,2,3 projectors, and 
\x is a singlet of SO(10)xSU(3) with VEV larger than Mgut- Equivalently one may try 
to directly generate the light Yukawa couplings as effective operators of the form 



A,; 



45^/45xY 



45- 



X 



LX 



Ix 



(4.2) 



12 3 

where the - ' are now antisextets of flavour and their VEV are the ^^1,2,3 projectors. 

All these effective operators can be realized using just renormalizable interactions 
introducing additional fields, that vary for the different cases. In general one has to add 
additional fermion multiplets beyond the 16' 16', that can be taken to be vector-like as 
well; they will be denoted as 16", 16", 16'", 16"', etc. We then allow couplings of all the 
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vector-like fermions together, but for simplicity we couple the light multiplet 16 just with 
16', 16', to realize the universal seesaw. The coupling matrix can be illustrated as: 





16' 


16[ 16" ... 


16^ 




10 








16" 




Mfuii 



and Mfuii will contain all the renormalizable couplings with SO(IO) or flavour fields such 
as 45, Xiji that all develop a VEV around the GUT scale. 

The condition to use the universal seesaw approximation is that the involved eigenval- 
ues of Mfj^ii are heavier than (45/j); once this is assured, the mass matrix Mheavy of the 
heavy states is effectively given by the projection of MJ^^ into the 3x3 subspace of the 



16^, 16;: 



A ~ (45h)Mj 



1 

heavy 



(45r)(1,0, 



M 



fall 







(4.3) 



where 1 stands for the identity in flavour space, 5ij. This means that to describe in the 
general case the resulting Yukawa couplings it is not necessary to diagonalize the entire 
coupling matrix, but it is sufficient to find the first 3x3 block of the inverse of Mf^u- 

In the following we describe some of the realizations of the Yukawa couplings ( [4.2| ) 
for the interesting cases n = 0, 1, 2, introducing the additional fermion multiplets and the 
appropriate flavon fields. 



4.1 Model I: first family unification, n = 

For n = the three terms appearing in A have negative powers, 0, -1, -2, of 45x/lx- 
These can be realized by introducing three more vector-like fermion multiplets that are 
flavour triplets, 16'/ 2 3 (with their SO(IO) conjugates 16'/ 2,3) and arranging the couplings 
as follows: 

16* 16^ 16'/* 16'2'* 16^'* 



16'/^' 
16'^^ 




where ^ 



1,2,3 



are three flavour (anti)sextets. The effective Aj of (4^) emerges here from 

1,2,3 



a double seesaw mechanism [^], and the VEVs of ^^j ' directly generates the rank-one 



flavour projectors: A,2,3 oc (C^'^''^)- Equivalently the Mheavy mass matrix of (|4.1|) and the 
0i,2,3 projectors can be seen to emerge from the first 3x3 block of the inverse of Mj^n. 
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Such realization is not unique, for example one may realize the same effective operators 
in a more symmetric fashion with less fields, by introducing just two vector-like triplets 
and antitriplets 16", 16-', 16"*, 16"* along with three sextets Xn , with couplings as follows: 



16^ 
16^ 
16; 
16"^ 



16* 


16^ 


16f 


16"* 




10 








Xi 








1 ij 
2X2 


ij 


Ix 




45x 







In this realization the VEV of the Xi,2.3 sextets gives the inverse flavour projectors that 
appear in Mhcavy^ i-e- Qi,2,3 oc (xi,2,3) in (|4.1|) . Then, the final Xf will still be given as a 
combination of -Pi, 2, 3 as from the inversion property ( 2.12| ). 

The coupling arrangments above, and the ones that will follow, are not the most general 
ones with the given fields: while some zeroes are motivated by SU(3) flavour symmetry 
or by SO(IO), some other choices require a further explanation (for example the choice 
between 45x or Ix)- This may be linked to the presence of some additional symmetry 
(e.g. Z3) that acts differently on the fields, or in the intertwining of gauge and flavour 
symmetries in a larger group. While these directions of investigation are very interesting, 
they go beyond the scope of the present work. 



4.2 Model II: second family unification, n = 1 

For the case n = 1, we describe first the realization using three triplets, that are equivalent 
to three sextets with rank-one VEV. One looks for a heavy mass Mhcavy of the form 



eavy 



CX 



(45x) 



(Ix) 



(4.4) 



We parametrize the three sextets Q12 3 as before via the tensor product of three triplets: 
Qn oc QuQln with the Qn the VEV of three flavour triplets scalar fields. Then the effective 
light Yukawa matrix can be realized by introducing four more vector-like fermion multiplets 
that are flavour singlets: 16", I62, I63 and 16"' (and their conjugates) and arranging the 
couplings as follows: 

16* 16^ 16'/ 16'2' 16^5' 16'" 
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By inverting in the vector-like sector the mass matrix Mj^n and projecting in the 16-, 16^ 
space one can see that the above form of Mheavy is reproduced, inverted. Therefore, via 
the universal seesaw, one obtains in the light sector the correct Yukawa matrices: 



A oc 



(45x> t t (Ix) t 



f = u,d, e, un 



(4.5) 



i.e. Pn oc PnPn^ with the reciprocal vectors defined again as Pn 



(as from the inversion property ( |2.12| )). 

It is interesting to note that the way to obtain direct and inverse powers of 45x/lx 
together is that the first two channels 16", I62 are realizing a "double" seesaw, while the 
third I63 a "triple" one. For we stress however that in the full seesaw the three rank-1 
contributions work together, to reproduce the three terms in A. 

One can also note that the scale of the ^1,2,3 triplets can be raised above the 45x scale 
and the double and triple seesaws still work as before. In the double seesaw this can be 
seen by comparing the determinant of the 2-3 block with that of the full matrix. A similar 
but more involved mechanism works for the triple seesaw. Indeed the next-to-smallest 
eigenvalues are, in the three cases, qi, q2/^x and 45x- In the third case this is valid 
when (/n > Ix > 45js:-^ Assuming then (7n > Ix > 45x ~ Mgut, one can obtain that 
no mass eigenvalue, apart from the light fermion masses, drops below MquTj so there will 
be no thresholds to take into account below the GUT scale. For this to happen, the scale 
of flavour breaking will be higher not only of the SO(IO) breaking scale, but also of the 
higher \x scale. 

Other realizations of this model can be built by using antisextets or, for example, two 
triplets qi^2 and one sextet ^3^ for this one must introduce three vector-like fermion fields: 
two flavour singlets 16/2 and one triplet 16"'* (and their SO(10)xSU(3) conjugates) and 
arranging the couplings as follows: 





16^ 


16^ 


16'/ 


16'^' 


16"'^ 


16^ 




10 








T6; 


45h 




qi 


92 


45x 


16'/ 




qi 


45x 






16^' 




92 




Ix 








X3 






Ix 



The Qi,2,3 projectors appearing in the effective operator will be given by Qi oc qiq\., 
Q2 oc 5292 ^'^'i O3 cx^ X3- Then the A, 2,3 projectors entering the light Yukawa matrices will 
be given again by the inversion formula ( p.l2| ). 

4.3 Model III: third family unification, n = 2 

The third case corresponds to the popular but troublesome unification of the third family 
quarks and leptons. It requires very large tan (3 that is currently disfavored, and also 
■^Considering that (10) < {45i{) < (45x>. 
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involves a nonlinearity also from the renormalization evolution of the large A;, Yukawa 



eigenvalue, as well as more evident corrections to the b mass |2C]. For these reasons we will 
not pursue the analysis of its phenomenological viability: we will just give an example of 
a possible realization of the effective operator in this last case. 

For n = 2 we look for a realization of the following effective operator: 



M, 



heavy 



CX 



(45 



(4.6) 



for which one is led to use a "triple" seesaw in all the three channels, by introducing for 
example six new flavour singlets multiplets 16'/ 2,3) 16i'2,3 (and their vector-like conjugates), 
with the following couplings: 





16* 


16^ 


16'/ 


16^' 


16'i 


16'/' 


16'2" 


16- 


16^ 




10 














T6; 






91 


92 


93 








16'; 




Qi 








45 






T6'2' 




Q2 










Ix 




T63 




Q3 












Ix 


I61 






45 






Ix 






I62 








45 






Ix 




I63 










Ix 






Ix 



5. Unification of determinants 



Case II above, with n = 1, is particularly interesting since it leads to an intriguing relation 
between Yukawa couplings. Looking at the determinant of the A/ matrices, it turns out to 
f actor ize as:* 

detXf = X^Dl, (5.1) 



with Dp = pi ■ {p2 Ap3). This expression is independent of ej, and this means that the 
determinants are unified among the four U, D, E, sectors. In other words, at GUT scale 
the eigenvalues will satisfy the following exact relation: 



Ae A^ At- — \d \s \b — A„ Ac At — Ajy^j ^vd2 ^t^DZ ■ 



(5.2) 



The symmetry structure encoded in the generation of the Yukawa matrices, for case II, 
ensures that this relation is satisfied at GUT scale independently of the model parameters 
and details such as the pattern of breaking in the flavour and gauge sectors. 

Obviously to rephrase (|5.2| ) as relations between fermion masses at the weak scale one 
has to take into account the RG running, and this reintroduces a dependence on the model 
details, the value of a^\z^ and mainly on the supersymmetry breaking scale. 



'One has in fact det(aiPi + a2-P2 + oizPz) ~ oi\a2azD^. 
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To test the consequences of ( |5.2| ) on the low energy masses we recall that, in the MSSM, 
fermion masses and mixing angles are defined as follows in terms of quantities at the GUT 
scale: 

XdV cos (3 TJudsRd 
XgV cos PlJudsRd 
XbV cos 13 VbRdBt 

= 0i2 B^^ , 



rriu 
rric - 
mt -- 



-- XuV sin (5 rjudsRuBf 
XcV sin P 7]cRuBt 
Xtv sin PrjtRuBf 



XeV cos (3 'nei_,Re 
f, - X^V cos P r]ef,Re 



XrV COS (3 rjr Re 



'12 1 Z 



712 



nslz 



hslz — Gi2 B. 



(5.3) 
(5.4) 



where: the factors Ru,d,e account for the running induced by the MSSM gauge sector, from 
Mgut to mt; the factor Bt accounts for the running induced by the large Xt; and r]e^^uds,b,c 
complete the QCD+QED running from mt down to 2 GeV for u, d, s or to the respective 
masses for b, c. 



From |21], assuming a-^lz = 0.118 and MgusY = "mt, we find 
r/„rf, = 1.74, Vc = 2.11, 776 = 1.52, ry, = 1.04, 77^^ = 1-06, 
Then, we calculate 



Ru = 3.45 



Rd = 3.36, 



Rp 



1.51 



R, 



1.39. 



(5.5) 



(5.6) 



The dependence on a^\z of all these coefficients affects mainly the lightest quarks via 
finds-, Vc (10%), all the others have a variation of the order of 1-3%. These uncertainties 
are however correlated. In the Ru,d factors there is also a stronger dependence on the 
super symmetry breaking scale. 

Finally, the factor Bt is a function of 
At, encoding the nonlinearity in the "top" 
renormalization due to its near Pendleton- 
Ross supersymmetric fixed point. We ob- 
tain Bt\x^.5 = 0.901, StU^^.s = 0.670. 
Near lower At the dependence on At is stiffer, 
and in power 6 it allows a quite precise 
determination of At from the experimen- 
tal value of mt- Indeed, one can estimate 
the factor Bt after determining At from mt, 
once one notices that mt is almost indepen- 
dent of tan P in its moderate range, i.e. sin P ~ 0.99 within 0.5%, and that At turns out to 
lie near 0.5, so that Bt is sufficiently dependent of Af. Taking into account the uncertainty 
in mt one finds At = 0.536-0.809 and Bt = 0.8734 ± 2.5%m,. 

At this point, identities ( |5.2| ) can be rewritten in the following form involving only low 
energy fermion masses: 




Figure 4: Bt as a function of At. 



rudms 



mem^mr RdVbV 



uds 



mb RivIfj,VT 

mem^mr RuVcVuds 
mjnt Rlril^Tlr 



Bt 



{tan PYBl 



12 



(5.7) 
(5.8) 
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and this implies the following predictions for rridms and m„: 

mams = 880 MeV^ ± 2%„, ± 4%„^ ± 25%„3 (5.9) 

m„ = f j 1.95 MeV ± 20%„, ± 10%„, ± 15%„3 . (5.10) 

Using a mean value of 19.5 for the known ratio ms/m^, the first equation gives the central 
values 

mrf = 6.7MeV, m^ = 131MeV, (5.11) 

that fit well the data, with a bit high, within la of the experimental range. This means 
that experimentally the D, E determinants are unified at GUT scale within 

The prediction for rriu carries a dependence on tan /?, that can be used to derive a 
prediction for it: 

tan /3 = 7.3 ± 1.3^„ ± 0.5^, ± 0.2^^ , (5.12) 
where we used the la range for = (3 it 1.5) MeV p2|. This prediction lies exactly in 



the favored range determined from the recent analysis in various MSSM scenarios 18|. 

An other prediction for tan /3 may be derived by combining the two above and assuming 
for mu/md and their experimental ranges. One gets: 

tan/3 - ( mcmA ^^^Ru f ybVuds V^^ 

= 10 ± 1.3^„/^, ± 1ms ± Im, ± 0.5„3 . (5.13) 

Of course this prediction tells what is the value of tan /3 once the model has been solved 
and rriu/md, mg and mt are determined. However we have found that the dependence 
on as 1 2 cancels almost perfectly in both the ratios involving the ii's and the r/'s, and 
enters only through the Bt dependence (< 4%). This is a remarkable fact, considering 
that for individual masses the uncertainty in as\z generates the dominant errors in their 
predictions. 

In addition, and equally remarkably, the even larger dependence on the supersymmetry 
breaking scale also cancels almost exactly in the ratio Ru/Rd, and enters just through the 
Bt factor, where it can be estimated as a small decrease of 3% when MsusY is raised from 
167 GeV to ITeV. 



6. Model analysis 

One may take the decomposition in rank-one projectors described until here as an ansatz, 
and verify that it can explain all fermion masses and mixing angles. To this aim we will 

®The same analysis carried out in the non-supersymmetric Standard Model gives an even better agree- 
ment of the D, E determinants (within 2%) at its 1-2 unification scale Mqu'p_]^2 — 10^'^ GeV. However, 
even ignoring the lack of simple unification of the three gauge coupling constants, the requirement of the 
U determinant to match with D, E requires an extension to models of the 2H[DM type, with the relative 
tan /3 parameter again in the moderate range. 
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write the mass matrices in terms of a number of independent parameters, and show how 
to determine them from the known masses, mixing angles, and CP violation phase. 

First of all let us parametrize the symmetric flavour projectors Pn in terms of the three 
triplets Pn, and choose a flavour basis to write the projectors as follows: 











fl a b\ 




/o 


0^ 




/o y 


A 




= X 




a c? ah 


V 


1 






000 








ah}? j 




c 


?) 


^00 J 



\ 



aef 
hef 



a 6/ 



1 + a^e 

c + ahe 



hef 

c + ahej 



(6.1) 



where we have rotated globally all the pn vectors in flavour space and rescaled e j , A to set 
pi = (1, a, 6), p2 = (0, 1, c), p3 = (0, 0, d). Eliminating also the irrelevant phases, we end up 
with the following 15 parameters domain: \,a,d £ M^, 6, c, £„, e^, Eg, e,y^ E C. Because of 
the SO(IO) relations ( p?]) two e may be eliminated and we remain with 11 real parameters. 
This form should account for all the Dirac couplings of charged fermions and neutrinos. 

For the Majorana mass matrix of neutrinos, in view of ( |2.29| ), we have a similar form, 
but involving e^,^ in place of and with a appearing in the d'^/ei,^ term in the 3-3 entry. 
The neutrino sector thus adds 4 '"ea/ parameters to the model. 

6.1 Leading order 

The yukawa sector can be analyzed quite completely by computing the leading order ex- 
pressions for eigenvalues and mixing angles. For charged fermions: 

\d^ 



X- 



A. ~ A. 



I + a^erfl 

II + a^eel 



Ac ~ A, 
Xs ^ X\l + a'^edl , 
X^ ~ A|l + a^eel , 



At- A 
Af,~ A 
A^ ~ A 



d^ 



012 ~ 



aed 



1 + a^ed 



O23 - ^ed 



7CP ^ arg 



Vd 



(6.2) 
(6.3) 
(6.4) 
(6.5) 



where rjd = {c + ahed)- Note also that since <C ed^e, the quark mixing angles are well 
approximated by the "down" sector. The analysis proceeds with the following steps: 

• First we note from the third generation that, to respect Xh/Xj- ~ 1, one should have 
I Eg I — \ed\ as already discussed. This relation holds within 10%, but is nevertheless 
important for the following first order calculation. 

• Then, if one notes that ignoring the corrections the second generation Yukawa are 
unified Ac,s,^ — A, on the other hand the RG invariant relation S12 — \JrndJrfis tells 



us that \c?td\ — 1) so that a should actually be quite large, \a\ 
later) and the corrections can not be neglected. 



\^d 



-1/2 



(~ 7, see 
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On the contrary, one sees that the required spUt of s,fj, Yukawa is produced by the 
factors Ad^e = |1 + a'^^d,e\ once the sign (the phases) of a^ee and a^e^ are nearly 
opposite, a^ee ^ —a^ed — e^-^\ so that \Ae\ ~ 2, \Ad\ ^ 0.5. For the first generation 
the spht goes (correctly!) in the opposite direction. Therefore, in this approach, the 
largeness of a, following from the Cabibbo angle, explains also the deviation from 
exactly family hierarchical masses of e, d, /i, s. This fact was observed in |^. 

From other known mass ratios one can then determine \€f,^d\ — 1/55, a ~ 7.5, je^l ~ 
1/600, ~ 0.5, tan/3 ~ 7.3. 

Further insight can be gained by looking at the other quark mixing angles. From: 

^13 ~ |6ei/(i2| =3-10-3, ^23 ~ I W*^'! =4-10-2, (6.6) 

where r]d = {c+abcd), one derives \b\ ~ 5, |c| ~ 0.5-^1.5. Therefore the consequences of 
^13 are similar to the effect of the Cabibbo mixing: while there we got that numerically 

— 1/2 —1/2 

a ~ , now we see that also & — is large, of the same order. This does not 

introduce significant correction to the third-generation eigenvalues but, below, will 
help neutrino to have a maximal atmospheric mixing. 



A phase can be determined from the CKM CP-violating phase 7cp 



.10 



7cp=^arg(^^i±^) =1.02, =^ arg(6) - arg(7?rf) = 1 ± 1 = { ° , (6.7) 

where we had from the first two steps that arg(l -|- a'^ed) — ±1, the sign depending 
on complex conjugation of e^. Fixing in this way the relative phase of b with rjd 
also fixes the modulus of c = r/^ — abcd- The two branches correspond to |c| ~ 0.2, 
|c| ~ 1.5. In the first branch, c and b are approximately orthogonal, in the second 
they are approximately opposite. 

• Two remaining free parameters can be taken to be |ee/erf|, argc, and one can be 
eliminated (using second order corrections to eigenvalues) from the precise ratio 
Xh/\r^[\ee/ed\+2Re{edcyd% 

Summing up, the leading order analysis for charged fermions shows a one parameter family 
of solutions, with two branches. The fiat direction can be taken as a combination of |ee/ed| 
with the common phase of 6, c. 

The precise leading-order solution for all parameters is not conclusive since next-to- 
leading corrections and contributions to CKM from the "up" sector are 10% in magnitude, 
larger than experimental errors; hence a numeric fit is necessary and will be described 
below. Nevertheless the model proves capable of accounting for all the charged fermion 
masses and mixings, giving a mechanism for linking the deviations A^^d with the Cabibbo 
angle 612- 



^°We assume that no other contribution to 7cp is introduced between GUT and our low energy scale. In 
such a case this formula should be modified accordingly. 
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We can now comment on the solution found from the charged fermions so far. Col- 
lecting the Pn vectors: 

pi = (1, 7, 5) , p2^ (0, 1, 0.2-1.5) , P3 - (0, 0, 0.7) , (6.8) 

we observe that the they tend to lie in the 2-3 plane. Exact planarity is broken to order 
1/8 ~ -y/ed, as can be seen from the first 1, with respect to the modulus of its vector 
\pi\ ~ 8.5. The reason for this is to be tracked in the magnitude of the quark mixing angles 
that require a and b to be large while c and d are of order one.^^ In terms of the reciprocal 



vectors (71,2,3 that appear in the model II realization in section [4.2| , quasi-planarity shows 
up as quasi- alignment. We have 

gi = (1,0,0), 52 = (7, 1,0), g3-(7,l,l). (6.9) 

We find this pattern of quasi-planar or quasi-aligned triplets a nice hint for the realiza- 
tion of in the fiavour sector of the theory, where three flavour triplets scalar fields acquire 
a VEV that generates the pn vectors. We leave this possibility for future model building; 



see e.g. [15| for realizations of disoriented VEVs from a potential. 

We observe also that if one may be concerned with the moduli of pn not being equal 
and of order 1, (they are in fact \pi\ ~ 8, \p2\ — 1.5, \ps\ ~ 0.7), they can be brought 
to be in the range 1-2 by rescaling all the e/ by a factor of 10. With this normalization 
the hierarchy of eigenvalues is due partially to the and partially to the pn vectors 
having hierarchic entries, while the normalization that we adopted above corresponds to 
an eigenvalue hierarchy that comes only from the e/'s. The choice may be varied depending 
on the theoretical realization in the fiavour sector, that should be the subject of a further 
study and goes beyond the scope of this work. 

6.2 Neutrino analytical 

From the light neutrino mass matrix (|]2|) and using the explicit parametrization given 
in (16.11), we can rewrite: 



V 

\aey/Ay\ \beu/A^\ 

\aeylAy\ 1 p 

\\heylAy\ p p^ + a 



with Ay = {1 + a^e^), p = {c + ahe^) / A^, a = {(fa + {h - acfel) / A^. 

In this matrix, we only have two free complex parameters and a, and the common 
phase of 5, c that was a flat direction of the charged fermions. However it is easy to see by 
homogeneity that this latter phase is irrelevant since it can be eliminated using the phase 
of a. Therefore, we have just two complex parameters and a, and also the neutrino 
sector is quite constrained. 



"^Aa a side note, the value found for S ~ 0.5 tells that the eigenvalues in the "up" sector as from 



line (6.2) do not strictly follow an exact hierarchy as in (1.3). A factor of 0.5 will decrease, with respect to 
precise hierarchy, either the top mass, which would be unacceptable, or the up quark mass, which is well 
tolerable. We will see also from the numerical solution that low niu is a prediction of this ansatz (see e.g. 



figure ^ . 
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Because of the large neutrino mixings and of the smah but important uncertainties in 
the parameters determined from the charged fermions, an exact analytical decomposition 
of this matrix is quite difficult. Nevertheless some useful observations can be made: 



First, since the ratio of the 1-1 with the 1-2, 1-3 entries is small ~ 1/7, from the 
known pattern of mixing angles it follows that this matrix predicts direct and non 
degenerate neutrinos: mi,i/mi,2 — 1/4 so that 111^2 /"mus — 1/5. 

Then, because the 1-2 entry is larger than the 1-3 one (1-2/1-3 ~ 1.5-2), then 0e3 
cannot be zero. 



• An estimate of the angles from the 2,3 sector can be given first in the form of a 
correlation between them, using the parameters found from the charged sector and 
after the requirement of maximal ^atm^ 



Oe3 _ my2 



1 + Re 



1 - Re^' 



With i?= - ~ ^i^-^ ~ ,_L^, (6.10) 
a 0;i2 A- 1.6-2.1' ^ ' 



where 6 = aigp. After imposing the correct neutrino hierarchy to solve for 6^/, we 
can estimate more explicitly: 



^soi ^ 30° 



1 -i^e 



iS 



0e3 - 6° 



1 + Re 



iS 



(6.11) 



Therefore, choosing e purely imaginary or just slightly negative, we see that the 
right ^soi — 32° is accommodated, with a prediction for Oes — 5°. 

The precise numbers are sensible to the precise values of the quark mixing angles 
and CP-phase, as well as to the neutrino hierarchy. In practice 6^3 can be brought 
reach the extrema of the range 0e3 — l°-7.5°, by stretching within la the angles and 
the neutrino hierarchy. This means that the model predicts a nonzero 9e3 inside the 
present 2a range (=0°-9.6°). With more precise experimental values of CKM angles 
and neutrino hierarchy, also the prediction for Ogs will become narrower. 

Finally, it can be seen that among the four parameters a flat direction emerges, along 
which the angles and hierarchy do not vary. It directly corresponds to the leptonic 
CP violation phase, that is thus left unpredicted by this ansatz, until angles and 
hierarchy will be known to better accuracy. Also one Major ana phase varies along 
this flat direction. 

Summing up, in the neutrino sector the model predicts hierarchic neutrinos with 
nonzero 9e3 in the present range, and there is a useful conspiracy from the charged fermions 
sector to allow the right neutrino pattern. A flat direction emerges as a combination of the 
four real parameters, and has to be added to the one found in the charged leptons sector. 
Therefore the whole model effectively takes advantage of just 13 of the 15 real parameters 
to reproduce the known data. 

For a precise test we performed a numerical analysis, where we fitted neutrinos together 
with charged fermions. The ^e3 mixing, that we have shown to lie automatically inside its 
allowed range, will not be used as a fit parameter but rather displayed as a model prediction. 
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6.3 Numeric fit and results 

The fit was performed using the 15 parameters (11 charged fermions, 4 neutrino) against 
20 tests for all known data ||2^, . The best fit results are shown in tables |^ and |2|, for 
models I and II. Each model has solutions in two branches, and model I shows difficulties 
in fitting the "up" sector (rn-t is too low, and rric too high). Model II instead can fit almost 
perfectly all data, albeit predicting a bit low rriu and a bit large nis as implied by the 
predictions described in section ^. 

The dependence on the supersymmetry breaking scale is the major source of theoretical 
uncertainty, and we have estimated its effect by performing the fit also with MsusY = 1 TeV 
and including mixed RG evolution. Raising the SUSY scale leads to better results in 
both models, but more appreciably in model II where all tension with experimental data 
disappears, and unification of determinants works perfectly. 

For model II we also illustrate the best fit Yukawa couplings in figure ^. 

A number of observations from the fit results in model II can be made: 

• Apart from and with lower MsusY that are predicted in the la range, all 
data are accomodated without any tension. Within the same total deviation one 
may choose even to adjust better either uis, rriu/md or rric and mt- 

• A flat direction and the two branches |c| ~ 0.2, 1.5 correctly emerge from the fit. In 
the second branch all complex phases of all parameters are almost aligned, hinting for 
a model with aligned efs and real parameters, with a reduction of the total number 
of real parameters to 10. 

• In the lepton sector, the mixing from terrestrial neutrinos is nonzero: Oez — 2°, 6°, 
respectively in the two branches. 

• Neutrino come out hierarchic and non degenerate niy ~ (0.0045, 0.01, 0.048) eV, and 
right-handed neutrino masses are of the order (10^ 10^ 10^^) GeV. 

• The leptogenesys asymmetry appears to have the correct sign ei ~ —3 • 10~^ at least 
in some branches, though (due to small RH neutrino masses) the value itself is too 
small for the usual thermal scenario, and it requires a speciflc assumption on the 
non-thermal production of the heavy neutrino initial abundances. 
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MODEL I 




exp 


1^' branch 


f = mt 
2"<^ branch 


-^SUSY 

I'** branch 


= ITeV 
2"<^ branch 



Fit 



mt 


167 ±5 


147 


148 


155 


155 


rrih 


4.25 ±0.15 


4.36 


4.38 


4.37 


4.37 


rric 


1.25 ± 0.10 


1.49 


1.47 


1.46 


1.45 




105 ± 25 


132 


119 


128 


116 


ms/md 


19.5 ±2.5 


21.1 


21 


21.8 


22.2 


EUipse 


23 ±2 


22.1 


22.6 


22.5 


22.9 


rriu/md 


0.5 ±0.2 


0.26 


0.27 


0.26 


0.25 


rriu+md 


8.5 ±2.5 


7.6 


8.0 


7.3 


6.5 




12.7 ±0.1 


12.75 


12.75 


12.81 


12.83 


102 023 


4.13 ±0.15 


4.11 


4.12 


4.11 


4.12 


103 013 


3.67 ±0.47 


3.82 


3.79 


3.82 


3.79 


sin2/3cp 


0.736 ±0.05 


0.735 


0.734 


0.730 


0.734 


7CP 


1.03 ±0.22 


0.93 


1.05 


0.93 


1.02 


^sol 


32 ±3 


31. 


32. 


31.0 


32. 


^atm 


45 ± 7 


45.3 


45. 


45.3 


45. 




0.18 ± 0.07 


0.105 


0.180 


0.171 


0.185 


Xdof 




0.68 


0.53 


0.49 


0.43 


Predictions 


nid (MeV) 


6±2 


6.0 


5.4 


5.9 


5.25 


Oe3 (°) 


- 9.6 


1.8 


6.6 


1.8 


7.2 


m„i (cV) 




0.0029 


0.0026 


0.0029 


0.0028 


(cV) 




0.052 


0.048 


0.052 


0.048 


TUee (eV) 




0.0007 


0.0009 


0.0007 


0.0009 


ei 




+2 10-12 


+6 10-12 


+210-12 


+2.710-12 


tan/J 




1.21 


1.23 


1.13 


1.14 



Table 1: Result of best fit against the known experimental ranges, for the two branches of model I, 
varying also the supersymmetry breaking scale between the "top" scale 167 GeV and 1 TeV. In bold 
are shown the quantities that show a la tension with experiment. Lepton masses are reproduced 

exactly and arc not displayed. Wc report also some predicted quantities, not used for fitting, but 
determined in the best fit solution. This model shows a evident difficulty in accommodating the 
"up" quark sector. 



-25- 



MODEL II 






AfsUSY = TO* 




= ITeV 




exp 


1^' branch 


2"<l branch 


branch 


2"<^ branch 


Fit 


mt 


167 ±5 


167.6 


169.0 


167.4 


167.8 




4.25 ±0.15 


4.40 


4.45 


4.32 


4.38 


rric 


1.25 ±0.10 


1.32 


1.31 


1.29 


1.29 


vis 


105 ±25 


134 


134 


127 


126 


rris/md 


19.5 ±2.5 


21.1 


21.4 


21.01 


21.0 


Ellipse 


23 ±2 


21.5 


22.1 


21.46 


22.2 


mu/nid 


0.5 ±0.2 


0.20 


0.25 


0.21 


0.32 


rriu+md 


8.5 ±2.5 


7.5 


8.0 


7.3 


7.9 


012 


12.7±0.1 


12.74 


12.77 


12.72 


12.75 




4.13 ±0.15 


4.12 


4.14 


4.11 


4.13 


10^ 013 


3.67 ±0.47 


3.77 


3.77 


3.77 


3.76 


sin2/3cp 


0.736 ±0.05 


0.730 


0.734 


0.73 


0.734 


7CP 


1.03 ± 0.22 


0.97 


1.07 


0.98 


1.06 


^sol 


32 ±3 


31.5 


32. 


31.4 


32. 


^atm 


45 ± 7 


45.1 


45. 


45.2 


45. 




0.18 ±0.07 


0.178 


0.185 


0.176 


0.185 


Xdof 




0.30 


0.27 


0.22 


0.14 


Predictions 


jud (McV) 


6±2 


6.4 


6.3 


6.1 


6.0 




- 9.6 


-1+2.2 




1.9 


5.9 


m^i (eV) 




0.0031 


0.0030 


0.0030 


0.0044 


rrivi (eV) 




0.050 


0.048 


0.050 


0.0481 


niee (eV) 




0.00068 


0.00091 


0.0007 


0.00085 


ei 




+310-1" 


-210-^ 


+210-1° 


-3.310-9 


tan/3 




7 1+0.5 
'•^-0.3 


7 n+0.5 
'•^-0.5 


7.7 


8.9 



Table 2: Result of best fit against the known experimental ranges, for the two branches of model II, 
varying the supersymmetry breaking scale between the "top" scale 167 GeV and 1 TeV. In bold are 
shown the quantities that show a la tension. For ^eS and tan (3 we also indicate the ranges obtained, 
by stretching these quantities up to an additional la variation. It can be noted that raising the 
supersymmetry scale to 1 TeV eliminates almost all tension, especially in the 2'^'^ branch. 
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tan 13 ^ 1 




: 10' 



: 10" 



: 10" 



Figure 5: Illustration of yukawa couplings at GUT scale for the three families of charged fermions, 
resulting from the best fit in model II for Msusy = 167 GeV, compared to the values calculated from 
experimental masses (as shown in figure Thin lines connect yukawa calculated from experimental 
masses, thick lines correspond to the best fit result. 
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7. Conclusions 

We presented a predictive model of SU(3) flavour symmetry in the context of SO(IO) 
supersymmetric GUT that can account for all the known fermion masses and mixings at 
less than la level, including neutrino data. The model has been analyzed with a leading- 
order analytical approach and with a complete numerical fit. 

In the quark sector, it predicts la low and, without raising the SUSY breaking 
scale, la large m^; all the other quantities are well accounted, including the CP phase. 
In the neutrino sector the model predicts direct and hierarchical neutrinos with non zero 
0e3 — 2°, 6° correlated to other quantities, but always within the present 99% c.l.. Also 
the correct sign of the leptogenesys asymmetry is an outcome of the numerical analysis, 
even though its magnitude is one or two orders of magnitude too low for the usual thermal 
leptogenesys scenario. 

The SO(IO) realization of this framework uses an "universal" seesaw mixing with super- 
heavy vector-like fermions to transfer their mass matrices to the light fermions. Together 
with adjoint Higgs fields, it has been shown to suppress the dominant part of D=5 pro- 
ton decay and supersymmetric flavour changing effects, while solving the doublet-triplet 
problem via the Missing VEV mechanism and allowing the correct mass generation. 

The flavour structure of the model is based on hierarchical combinations of three rank- 
one projectors in flavour space, built with three flavour triplets or sextets, that break the 
SU(3) symmetry. As opposed to traditional unification schemes, this model does not unify 
some Yukawa matrices at GUT, but realizes in an interesting way the exact unification of 
their determinants. 

A quite robust prediction for tan /? ~ 8 in the moderate, favored, range is also conse- 
quence of this kind of unification. 

The solution emerging from the analysis points to an interesting configuration of quasi- 
aligned flavour triplets, that suggests a precise pattern of SU(3) breaking. We leave this 
analysis for further work. Finally, the intertwining between flavour and SO(IO) Higgs fields 
also suggests that a similar mechanism in the context of higher unification groups may lead 
to a framework as successful as this one. 
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